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Improved flavor symmetry in Kogut-Susskind fermion actions 
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We study improved Kogut-Susskind fermion actions focusing on flavor symmetry restoration. Several variants 
of fat actions suitable for dynamical simulations are considered, including an action with no tree level 0(a 2 ) 
errors. The spectrum of all the pions is computed and used as a measure of flavor symmetry violation. Finally, 
the Naik term is introduced to restore rotational symmetry. 



1. INTRODUCTION 

Dynamical fermion simulations are still one of 
the challenging problems in lattice QCD. Re- 
cently, a lot of work has been devoted to devel- 
oping improved lattice fermion actions. In this 
paper, we report our findings on improved Kogut- 
Susskind fermions. In particular, we focus on the 
issue of flavor symmetry, and the construction of 
an action suitable for dynamical simulations. 

Kogut-Susskind (KS) fermions have 0(a 2 ) lat- 
tice artifacts, in contrast to the Wilson fermions 
which have 0(a). They represent four degener- 
ate quark flavors with a £7(4) x £7(4) symmetry 
in the continuum. This symmetry is broken down 
to a £7(1) x £7(1) on the lattice. This remnant 
symmetry protects against additive mass renor- 
malizations. The C/(4) x 11(4) symmetry break- 
ing manifests itself in mass splittings in the pion 
spectrum. There is only one Goldstone pion asso- 
ciated with the exact U(l) chiral symmetry, while 
the rest of the fifteen pions E| are non-Goldstone 
massive particles. 

The first step towards improving the KS ac- 
tion was made by Naik. He introduced the so 
called "Naik term" (— gA^), which eliminates a 
tree level 0(a 2 ) artifact, improving the Lorenz 
symmetry (dispersion relation). The remaining 
order 0(a 2 ) artifacts are the £7(4) x £7(4) break- 
ing terms. As Lepage pointed out, these terms 
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are due to couplings to high momentum gluons, 
which cause scattering of lattice quarks among 
the corners of the Brillouin zone. These scat- 
terings correspond to spin-flavor transformations 
of the physical quark field, causing the breaking 
of the symmetry. Thus, the way to reduce the 
symmetry breaking is to smear P-pL|| the quark- 
gluon vertex, which suppresses the couplings to 
high momentum gluons. Recently, it was shown 
by Lepage B that the one gluon exchange in- 
teraction at high momentum is identical to 4- 
fermi £7(4) x £7(4) breaking terms. Consequently, 
these effects can be canceled by either appropri- 
ate smearing of the quark-gluon vertex, so that 
high momentum gluons decouple, or by including 
in the action appropriate 4-fermi terms. 



2. ACTION CONSTRUCTION 

Since the £7(4) x £7(4) symmetry breaking oc- 
curs due to coupling to high momentum gluons, 
we need to modify the quark-gluon vertex in or- 
der to decouple them. The coupling of quarks to 
a gluon with any of its transverse momen- 
tum components p v = n/a causes quark scat- 
tering among the corners of the Brillouin zone. 
Couplings to gluons with longitudinal momentum 
p^ = -k / a vanish due to cancellations between the 
forward and the backward part of the Dirac oper- 
ator, consequently we do not need to worry about 
them. There arc 3 independent couplings to high 
transverse momentum gluons (Vi , Vi , V3) . We can 
eliminate them by modifying the gauge link used 
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in the lattice derivative as following: 

V 

P a 

S$(x) = U^xp^x + y)Ut{x + A) 

S$ p (x) = U v (x)S$(x + v)Ut(x + A) 

Sft^x) = U v (x)Sj$,(x + 0)U}(x + A) (2) 

It is not difficult to see that the above set of 
terms is the simplest needed to eliminate com- 
pletely couplings to high momentum gluons. 

In the weak coupling limit, the couplings 
(^11^2, V3) to the gauge field with one, two or 
three of the transverse momentum components 
±7r/a can be expressed as functions of the staple 
couplings W3 , , W7 , and the single link coupling 
ci: 

Vi = ci + 21U3 — 8w$ — 48ui7 
V2 = c\ — 2w3 — 81U5 + 4&W7 
V 3 = a- 6iu 3 + 24w 5 - 48w 7 (3) 

The overall normalization condition 

ci + 6iu 3 + 24w 5 + 48w 7 = 1 (4) 

is used to ensure that the total coupling to the 
nearest neighbor in the free field limit is one. The 
solution to the equations V\ = 0, V2 = 0, = 
is ci = 2u>3 = 8u>5 = 48w7 = 1/8. This set of pa- 
rameters defines our "Fat7" action. The "Fat5" 
action is constructed by using paths up to length 
5, minimizing the maximum of the couplings \V\. 
The "Fat5" couplings are c\ = 2w 3 = 8u>5 = 1/7, 
Wj = 0, which give \V\\ = V2 = | V3 1 = 1/7. 
Adding the Naik term to the "Fat 7" action results 
in an action with no tree level 0(a 2 ) errors and 
improved flavor symmetry. Unfortunately, the ex- 
tra couplings introduced to fix the high momen- 
tum behavior introduce a low momentum artifact 
~ A^kla 2 . In order to correct this new artifact, 
one has to introduce a new term discovered by 
Lepage H . This new term is a double staple 

PHHx) = U u (x)S$(x + v)Ut(x + A) (5) 



with weight vol- Cancellation of the new artifact 
occurs, if W3 — > W3 — 6wl and wl = —1/16. Note 
that because the new term couples to a distance 2 
in the v direction, at momentum p u — 7r / 'a it acts 
as a single link term. For this reason, the above 
change of couplings preserves the normalization 
condition and keeps the couplings Vi , V%, and V3 
zero. Adding the Naik term, (— ^A^) we get an 
action which we call "Asq", that has errors of 
0(a 4 , g 2 a 2 ). All the above actions can be tadpole 
improved. For the tadpole improved actions, we 
append the suffix "tad" to the names above. 

Motivated by the success of RG inspired APE 
smeared actions in the work of DeGrand, Hasen- 
fratz and collaborators || , we also decided to ap- 
proximate APE smeared actions as a set of paths 
with appropriate weights. This is necessary, if 
we want to use such actions in dynamical simula- 
tions. In |Q] we show how this approximation can 
be done expanding in powers of the APE smear- 
ing parameter a. For one APE smearing with 
a = .25, and to leading order in a we get the ac- 
tion we call "OFUN25" (once fattened unitarized 
a = .25). The fat link for this action is: 

- U M + f £ - U.S^l V U») . (6) 

3. SIMULATIONS AND RESULTS 

Our action tests are focused on the flavor sym- 
metry restoration, since we have established in 
our previous studies Q that the Naik term im- 
proves Lorenz symmetry adequately for the lat- 
tice spacings used in current simulations. Fur- 
thermore, the Naik term has little effect on flavor 
symmetry. Thus, in most of the smearing vari- 
ants we use, we do not include the Naik term. 
We include the Naik term into "Asq" and "Asq- 
tad" , since these actions are the candidates for an 
action with both improved flavor and rotational 
symmetry. The flavor symmetry breaking is es- 
timated by the splittings in the pion spectrum. 
Spectroscopy is done on the same set of lattices 
(sample sizes range from 48 to 60 ) for all the 
actions tested. We use 12 3 x 32 (3 = 7.3 lattices 
with two dynamical quark masses m = .020 and 
m = .040, and 16 3 x 48 (3 = 7.5 lattices with dy- 
namical quark masses m — .015 and m = .030. 
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Figure 1. Interpolated pion masses for all the actions tested. The highest level is the flavor singlet pion 
while the lower is the Goldstone pion. The first doublet is the local non-Goldstone pion (flavor structure 
7075) and the 7^75 pion (right). The second is the 7o7i(left) and 7^ pion. The third is the 70 (left) and 
7i Pion. 



The gauge action is 1-loop Symanzik improved 
and the sea quark action is Kogut-Susskind with 
a staple and a Naik term . We also use a set 
of 32 3 x 64 quenched lattices at = 6. 15 with Wil- 
son gauge action. The pion spectra are computed 
for the two different masses in the case of the dy- 
namical lattices and for m = .010 and m = 0.020 
in the case of the quenched lattices. We then in- 
terpolate the results to m w /m p — .55. The inter- 
polated spectra are shown in Figure [jj. In this fig- 
ure, we also show for comparison the spectrum of 
the standard KS action (OL) and the spectrum of 
the action used for the dynamical quarks (OFN). 
Our results in Figure [l] show that there is a sig- 
nificant improvement of the flavor symmetry as 
we increase the suppression of the high momen- 
tum gluon couplings. The highest improvement 
is observed for the "Fat7tad" and "Asqtad" ac- 
tions where couplings are set to zero at tree level. 
For all the actions the tadpole variant is better 
than the non-tadpole. The approximately reuni- 
tarized action "OFUN25" is slightly better than 
the "OFN" . In all cases the spectrum exhibits the 
degeneracies suggested in M. For a detailed er- 
ror analysis that better justifies our claims see our 
recent publication M. 



Although the differences in flavor symmetry 
improvement among the various highly smeared 
actions seem small, we favor "Asqtad" on theo- 
retical grounds (0(a 4 , g 2 a 2 ) errors), and on our 
estimation that it only costs a factor of 2-3 for 
dynamical simulations with light quarks. We be- 
lieve that such an action will be useful in realistic 
projects with dynamical quarks. 
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